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ABSTRACT 

Utilizing 21 new Chandra observations as well as archival Chandra, ROSAT, and XMM-Newton data, we 
study the X-ray properties of a representative sample of 59 of the most optically luminous quasars in the 
Universe (M, w —29.3 to —30.2) spanning a redshift range of z ~ 1.5—4.5. Our full sample consists of 32 
quasars from the Sloan Digital Sky Survey (SDSS) Data Release 3 (DR3) quasar catalog, two additional objects 
in the DR3 area that were missed by the SDSS selection criteria, and 25 comparably luminous quasars at z>4. 
This is the largest X-ray study of such luminous quasars to date. By jointly fitting the X-ray spectra of our 
sample quasars, excluding radio-loud and broad absorption line (BAL) objects, we find a mean X-ray power- 
law photon index of F = 1 .92+^ and constrain any neutral intrinsic absorbing material to have a mean column 
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density of Nh ~ 2 x 10 cm . We find, consistent with other studies, that F does not change with redshift, 
and we constrain the amount of allowed F evolution for the most-luminous quasars. Our sample, excluding 
radio-loud and BAL quasars, has a mean X-ray-to-optical spectral slope of a ox = —1.80 ±0.02, as well as 
no significant evolution of Oo X with redshift. We also comment upon the X-ray properties of a number of 
notable quasars, including an X-ray weak quasar with several strong narrow absorption-line systems, a mildly 
radio-loud BAL quasar, and a well-studied gravitationally lensed quasar. 

Subject headings: Galaxies: Active: Nuclei — Galaxies: Active: Optical/UV/X-ray — Galaxies: Active: 
Evolution — Methods: Statistical 



1. INTRODUCTION 

The most optically luminous known quasars (M, < —29) 
serve as valuable astrophysical probes of extreme accretion 
conditions and the distant universe. These objects have been 
found to date at z ~ 1.5-4.5; their resulting i magnitudes of 
ps 15-18 and relatively bright multiwavelength fluxes allow 
them to be studied effectively with a variety of facilities across 
the electromagnetic spectrum. Even if they are radiating near 
the Eddington limit with L/Ledd ~ 1> their energy outputs re- 
quire ps 10 9 — 10 10 M nuclear black holes and thus they are 
presumably associated with the most-massive galaxies; to- 
day many of these objects have likely evolved into super- 
giant ellipticals found in the cores of rich clusters. As the 
most-luminous, non-transient emitters at high redshift, these 
quasars have been useful in cosmological studies including 
measuring absorption lines from intervening line-of-sight ma- 
terial (e.g., Rauch 1998; Wolfe et al. 2005; and references 
therein), assessing the Cold Dark Matter cosmogony (e.g., 
Efstathiou & Rees 1988; Turner 1991; Springel et al. 2005), 
and constraining the accretion history of the Universe (e.g., 
Croom et al. 2004; Richards et al. 2006). 

The Sloan Digital Sky Survey (SDSS; York et al. 2000) is 
now providing the most-complete selection of highly lumi- 
nous quasars to date (e.g., Schneider et al. 2005, hereafter 
S05). A large fraction of these objects, about 2/3, lack pointed 
or serendipitous X-ray detections (aside from at z > 4, where 
pointed X-ray observations have detected a large fraction; 
e.g., Vignali et al. 2003, 2005). Accordingly, we have started a 
project aimed at improving understanding of the X-ray prop- 
erties of the most-luminous known quasars over as broad a 
redshift range as possible. The X-ray emission from quasars 
probes the innermost regions of their accretion-disk coronae 
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where any changes in the mode of accretion might be most 
evident, and X-ray spectroscopy provides constraints on in- 
trinsic and intervening absorption. 

Recent studies of the X-ray spectra (e.g., Page et al. 2005; 
Shemmer et al. 2005a, 2006a; Vignali et al. 2005) and X- 
ray-to-optical spectral energy distributions (SEDs; e.g., Strat- 
eva et al. 2005; Steffen et al. 2006, hereafter S06) of quasars 
have generally shown no clear changes with redshift, although 
some exceptions have been found (e.g., Grupe et al. 2006; 
Kelly et al. 2007) and at lower luminosities X-ray spectral 
evolution may be observed (Dai et al. 2004). These results in- 
dicate that the inner regions of quasars are largely insensitive 
to the enormous changes in large-scale cosmic environment 
occurring over the history of the Universe. There is evidence, 
however, that the photon index (r) of the X-ray power-law 
spectrum increases as L/Lgdd increases (e.g., Shemmer et al. 
2006b) and that the X-ray-to-optical flux ratio (a ox ; Tanan- 
baum et al. 1979) drops with increasing luminosity. Studies 
of the luminosity and redshift dependence of F and a ox benefit 
from the widest possible sampling of the luminosity-redshift 
plane; such wide coverage is needed to break the luminosity- 
redshift degeneracy invariably present in flux-limited sam- 
ples. By systematically studying the most-luminous quasars 
over the full redshift range where they exist in the Universe, 
zsi 1.5-4.5, it is possible to populate one important region 
of this plane, complementing efforts to fill other regions of 
luminosity-redshift space (e.g., S06). Systematic X-ray mea- 
surements of the most-luminous quasars also serve to broaden 
the well-sampled luminosity range available for study and 
thereby minimize the possibility of confusion by spurious cor- 
relations (e.g., Yuan et al. 1998). 

In this paper we study, using a combination of new Chandra 
"snapshot" observations as well as archival Chandra, ROSAT, 
and XMM-Newton data, the basic X-ray properties of 32 of the 
33 most-luminous quasars in the SDSS Data Release 3 (DR3) 
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quasar catalog (S05; see §2.1.1 for a discussion of the one 
quasar that is not included in our study). All 32 of the quasars 
in our SDSS sample have X-ray detections. We also include 
two comparably luminous quasars missed by the SDSS se- 
lection and an additional 25 comparably luminous non-DR3 
quasars at z>4. We use our results to strengthen constraints 
upon the X-ray spectral and X-ray-to-optical SED properties 
of the most-luminous quasars, via a combination of single- 
object and multiple-object analyses. 

We detail the general properties of our sample in §2, as 
well as the X-ray observations and data reduction. X-ray, 
optical, and radio properties are presented in §3, and opti- 
cal spectra and notes on exceptional objects appear in §4. 
Data analysis and results are given in §5, and a summary 
of our findings is given in §6. We adopt a cosmology with 
H = 7010ns- 1 Mpc- 1 , i2 M =0.3, and£2 A = 0.7. 

2. SAMPLE AND X-RAY DATA 
2.1. Sample Selection and Properties 
2.1.1. SDSS DR3 Quasars 

About half of our sample of highly luminous quasars has 
been drawn from the SDSS DR3 quasar catalog (S05). The 
SDSS, an optical imaging and spectroscopic survey that aims 
to cover about one-quarter of the entire sky, targets active 
galaxies for follow-up spectroscopy primarily based upon 
their ugriz (Fukugita et al. 1996) colors and magnitudes (e.g., 
Richards et al. 2002). Active-galaxy candidates at z<3 are 
spectroscopically targeted if their i magnitudes are 15-19.1; 
high-redshift candidates are targeted if i = 15-20.2 (the limit 
at i = 15 is imposed to avoid saturation and fiber cross-talk 
problems in the SDSS spectroscopic observations). The DR3 
quasar catalog has been constructed from SDSS spectroscopic 
observations over a solid angle of 4188 deg 2 (about 10% of 
the sky). Given the large areal coverage, this catalog should 
contain representative members of the population of the most 
optically luminous quasars in the Universe; i.e., other surveys 
are unlikely to find a population of quasars significantly more 
luminous than those studied here. 2 About 60% of the most 
optically luminous quasars in the SDSS DR3 quasar catalog 
had been discovered in earlier surveys, such as the Hamburg 
Quasar Survey (e.g., Hagen et al. 1999), the Second Byurakan 
Survey (e.g., Stepanian et al. 2001), the University of Michi- 
gan Survey (e.g., MacAlpine & Lewis 1978), and the Palomar 
Digital Sky Survey (DPOSS; e.g., Djorgovski et al. 1998). 

We sorted the DR3 quasar catalog upon M, and consid- 
ered the 33 most-luminous quasars in the catalog for Chan- 
dra targeting (see Figure [Tjl. The number 33 was chosen 
based upon practical X-ray observing-time considerations, 
and this sample size is large enough to provide statistically 
meaningful results. Of the 33 most-luminous quasars, 1 1 al- 
ready had detections in archival X-ray data and were not tar- 
geted; these archival data have been utilized in our study. 
The remaining 22 quasars were proposed via the Chandra 
Cycle 7 Guaranteed Time Observing program, and 21 of 
them were awarded observing time. One of our targets, 
SDSS J10071 1.81+053208.9, was awarded to another Chan- 
dra observer (S.F. Anderson) as part of a program studying 
bright and extreme Broad Absorption Line (BAL) quasars. 
We do not consider the omission of SDSS J1007+0532 from 
our sample to be statistically problematic. In fact, owing to 

2 The recently released SDSS DR5 quasar catalog (Schneider et al. 2007), 
covering 5740 deg 2 , further supports this assertion. 




FIG. 1. — Absolute i-band magnitude vs. redshift for our SDSS sam- 
ple compared with the SDSS DR3 quasar catalog. Our SDSS sample 
of 32 quasars includes both archival (open symbols) and targeted (filled 
symbols) sources with Chandra (circles), XMM-Newton (triangles), and 
ROSAT (squares) observations. The two additional sources that were missed 
by the SDSS (APM 08279+5255 and HS 1603+3820; see §2.1.2) are 
shown as stars. The gravitationally lensed quasars APM 08279+5255, 
SDSS J0145-0945, and SDSS J0813+2545 have been de-amplified to their 
true Mj magnitudes and have bold symbols for clarity; all fail our cutoff at 
Mj < -29.28, which is shown as a dashed line. SDSS J1007+0532, which 
has been targeted by S. F. Anderson and thus needed to be removed from our 
sample (see §2.1.1), is shown as an open diamond. Small dots represent the 
ks 46,000 quasars in the DR3 catalog. 

its BAL-quasar nature, this object would need to be removed 
from most of our analyses of a ox , T, and other properties in 
any case. Our SDSS sample thus includes 32 quasars with M,- 
values of —29.28 to —30.24, all of which have sensitive X-ray 
coverage; we adopt M,- = —29.28 as a practical minimum lu- 
minosity for our sample. These 32 quasars span essentially 
the entire range of redshift (z«l .5^4.5) over which such lu- 
minous objects are known, although the source statistics at 
z > 4 are limited. 

2.1.2. Incompleteness and Complementary z^4 Quasars 

For the highly luminous and optically bright broad-line 
quasars under consideration here, the SDSS is not expected 
to suffer from substantial incompleteness biases. However, 
some incompleteness is expected due to the SDSS spectro- 
scopic limit of i = 15 and the fact that, at 2.5-2.9 and 
Z« 3.4-3.6, the SDSS colors of quasars intersect the stel- 
lar locus (e.g., Richards et al. 2002, 2006; S05). Further- 
more, about 5% of quasars are expected to be missed by 
the SDSS, largely due to image defects and source blend- 
ing (Vanden Berk et al. 2005). We have searched for missed 
quasars more luminous than M, = —29.28 in the area cov- 
ered by the DR3 quasar catalog using the NASA Extragalac- 
tic Database 3 (NED) combined with accurate SDSS photom- 
etry (since the photometry in NED is not uniform and can 
contain significant errors). We have found only two missed 
quasars: APM 08279+5255 at z = 3.91 (e.g., Irwin et al. 

1998) and HS 1603+3820 at z = 2.51 (e.g., Dobrzycki et al. 

1999) . 4 APM 08279+5255 is a gravitationally lensed BAL 
quasar with i = 14.9 that slightly violated the SDSS spec- 

3 See http://nedwww.ipac.caltech.edu/ 

4 Our NED search also uncovered the object 
QUEST J150724.0-020212.8 in the SDSS area, which is identified as 
an exceptionally luminous z = 1.09 quasar by Rengstorf et al. (2004). 
However, that identification relied on the assumption that the single feature 
seen in a spectrum was the Mg II line. Given this tenuous identification we 
obtained a spectrum of this source using the Low Resolution Spectrograph 
(LRS) on the Hobby-Eberly Telescope (HET; Ramsey et al. 1998). Our 
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troscopic limit. HS 1603+3820 has a very rich spectrum 
of associated heavy element absorbers (e.g., Misawa et al. 
2005); it is near the redshift range where the SDSS col- 
ors of quasars intersect the stellar locus and was missed by 
the SDSS quasar color-selection algorithm by just 0.02 mag 
in u (G. T. Richards 2006, private communication). For- 
tunately, both APM 08279+5255 (Chartas et al. 2002) and 
HS 1603+3820 (PI: Dobrzycki) have X-ray detections in 
archival pointed Chandra observations, and we include these 
quasars in our analyses below as appropriate. After these two 
additions, we expect < 15% incompleteness forM, < —29.28 
broad-line quasars at z ~ 1.5—4.5 in the area covered by the 
SDSS DR3 quasar catalog. The observation log of our sam- 
ple, including these two additional sources, appears in Ta- 
ble 1 ; this 34 object sample will hereafter be referred to as 
the "core" sample. 

In some analyses below, we will complement our core 
sample with 25 additional comparably luminous (M,<— 29) 
quasars at z>4 that lie outside the area covered by the SDSS 
DR3 quasar catalog. The X-ray properties of these quasars 
have been studied by Vignali et al. (2003, 2005) and Shemmer 
et al. (2005a) and appear in Table 2. These quasars were dis- 
covered by the DPOSS and other large-area surveys, and their 
basic optical properties are consistent with those of our SDSS 
quasars. This complementary sample is required to improve 
our statistical coverage at z > 4, which would otherwise be 
poor. Since the X-ray and optical properties of these objects 
were investigated as a whole in S06, this subsample will be re- 
ferred to as the complementary high-luminosity z> 4 quasars 
from S06. We have not added highly luminous quasars out- 
side the SDSS DR3 area at z < 4 from other surveys, since 
such quasars generally do not have systematic sensitive X-ray 
coverage. 

2.1.3. Radio-Loud, Broad Absorption Line, and Lensed Quasars 

Radio-loud quasars (RLQs) are known to have jet-linked 
X-ray emission components that generally lead to higher 
X-ray-to-optical luminosity ratios than those of radio-quiet 
quasars (RQQs; e.g., Worrall et al. 1987). Given this find- 
ing, it is important to consider RLQs and RQQs separately in 
statistical analyses of quasar X-ray properties. We quantify 
radio loudness using the radio-loudness parameter, R, defined 
as R = /5GHz//44ooA (Kellermann et al. 1989). We classify any 
quasar with R > 10 as radio-loud; details of our R calculations 
are given in §3. Two of the quasars in our 34 object core sam- 
ple, SDSS J0844+0503 and SDSS J1733+5400, are RLQs. 
Both objects are only mildly radio loud (R w 19 and R ps 10, 
respectively) and were observed as part of our Chandra Cy- 
cle 7 observations. Given the radio coverage of all of our 
sources (including tight upper limits of R < 0.5-4 on many of 
them), we do not expect there to be any RLQs that have not 
been identified as such. Although it is possible that there may 
still be some jet-linked X-ray contribution for sources with 
R < 10, we use equation (2) of Worrall et al. (1987) to esti- 
mate a limit of <2% on the fraction of X-ray emission from 
our RQQs that is jet related. 

Owing to intrinsic X-ray absorption, BAL quasars also re- 
quire special consideration in statistical analyses of quasar 
X-ray properties (e.g., Gallagher et al. 2006 and references 
therein). Using a catalog of BAL quasars from DR3 com- 
piled by Trump et al. (2006), we found six potential BAL 



quasars in our SDSS sample. The three we admit as BAL 
quasars are SDSS J0844+0503, SDSS J1525+5136, and 
SDSS J2313+0034. All three have positive "balnicity" in- 
dices (see Trump et al. 2006); note that SDSS J0844+0503 
is a mildly radio-loud BAL quasar. The other three potential 
BAL quasars, SDSS J1001+5454, SDSS J1407+6454, and 
SDSS J 1426+6025, have UV absorption but do not formally 
satisfy the positive balnicity criterion, so we do not remove 
them from our analyses below. The removal of these three 
sources from our main analyses does not significantly affect 
any of the results. Three of these quasars (SDSS J0844+0503, 
SDSS J1001+5454, and SDSS J1407+6454) were observed 
as part of our Chandra Cycle 7 observations; the remaining 
ones have archival X-ray coverage. Given that the redshifts 
of our quasars are sufficiently high (with the possible excep- 
tion of SDSS J0813+2545) so that the definitive C IV BAL 
transition lies within the spectral range covered by the SDSS 
spectra (this requires z> 1 .5), we do not expect there to be any 
unidentified BAL quasars within our sample. Note that the 
complementary high-luminosity z>4 quasars from S06 have 
been chosen to be radio-quiet, non-BAL quasars. 

Three quasars with archival X-ray observations described 
in this paper, SDSS J0145-0945, SDSS J0813+2545, 
and APM 08279+5255, are gravitationally lensed. 
SDSS J0145-0945 and APM 08279+5255 have flux- 
amplification factors of w 3 and « 100, respectively (e.g., 
Surdej et al. 1988; Egami et al. 2000; Lehar et al. 2000; E. O. 
Ofek 2006, private communication), and we have calculated 
the flux-amplification factor of SDSS J0813+2545 to be 
« 6 (based on the V-magnitudes taken from the discovery 
paper of Reimers et al. 2002). After correcting for flux 
amplification due to lensing, none of these quasars satisfies 
our Mj = —29.28 cutoff. Therefore, we present the basic 
X-ray properties of these three quasars below but exclude 
them from most of our statistical analyses. Since our sample 
consists of sources at the top of the quasar luminosity func- 
tion (resulting in a strong magnification bias), the expected 
fraction of lensed quasars is of the order of a few percent 
(e.g., Turner, Ostriker, & Gott 1984; E. O. Ofek 2007, private 
communication); this is consistent with the ~ 10% fraction of 
lensed quasars that we find for our highly luminous sample. 
We do not expect unresolved (< 1") gravitational lenses to be 
affecting our results materially, and we constrain extended 
X-ray emission in §2.2. 

2.2. Chandra Observations and Data Reduction 

Our 2 1 Chandra Cycle 7 targets (see §2.1.1) were observed 
using the Advanced CCD Imaging Spectrometer (ACIS; 
Garmire et al. 2003) with the aimpoint on the S3 CCD. 
The requested "snapshot" exposure for each target was 4 ks. 
All targets were placed near the aimpoint; with the ex- 
ceptions of SDSS J1350+5716, SDSS J1421+4633, and 
SDSS J1521+5202, all were strongly detected with w 10-150 
counts from 0.5-8 keV (details on SDSS J1521+5202 appear 
in §4). Faint mode was used for the event telemetry format, 
and all observations were free from strong background flares. 

Data analysis was carried out using standard CIAO 5 V3.2 
routines, and only events with ASCA grades 0, 2, 3, 4, and 6 
were used. WAVDETECT (Freeman et al. 2002) was used for 
source detection with wavelet scales of 1, \/2, 2, 2\/2, and 4 
pixels. We adopted a WAVDETECT false-positive probability 



spectrum does not show the emission feature seen by Rengstorf et al., and 
we identify this object as a Galactic star and not a quasar. 



5 Chandra I nteractive Analysis of Observations. 
http://asc.harvai-d.edu/ciao/ 
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TABLE 1 

X-ray Observation Log of the Core Sample 



X-ray Chandra Exp. Time a Alternate 



Object (SDSS J) 


z 


Obs. Date 


Cycle No. 


(ks) 


Designation 


Notes b 


012156.04+144823.9 


2.87 


2006 Jun 21 


7 


3.9 


HS 0119+1432 


PC 


014516.59-094517.3 


2.73 


2004 Aug 23 


5 


13.3 


UM673 


AC, Lensed 


020950.71-000506.4 


2.85 


2005 Dec 02 


7 


2.2 


UM 402 


PC 


073502.31+265911.4 


1.97 


2005 Dec 05 


7 


4.0 




PC 


075054.64+425219.2 


1.90 


2006 Sep 05 


7 


4.0 


HS 0747+4259 


PC 


080342.04+302254.6 


2.03 


2006 Sep 06 


7 


4.0 


HS 0800+3031 


PC 


081331.28+254503.0 


1.51 


2002 Jan 30 


3 


5.0 


HS 0810+2554 


AC, Lensed 


084401.95+050357.9 


3.35 


2006 Feb 18 


7 


3.8 




PC, BAL, RLQ 


090033.49+421546.8 


3.29 


2006 Feb 09 


7 


3.9 


HS 0857+4227 


PC 


094202.04+042244.5 


3.28 


2006 Feb 08 


7 


4.1 




PC 


095014.05+580136.5 


3.96 


2006 Jan 29 


7 


3.9 


PSS J0950+5801 


PC 


100129.64+545438.0 


1.76 


2006 Jan 29 


7 


4.0 


Mrk 132 


PC 


101447.18+430030.1 


3.13 


2006 Jun 14 


7 


4.1 


HS 1011+4315 


PC 


110610.73+640009.6 


2.20 


2006 Jul 16 


7 


3.7 


HS 1103+6416 


PC 


111038.64+483115.6 


2.96 


2002 Apr 25 




14.1 


Ql 107+487 


AX, Page et al. (2005) 


121930.77+494052.3 


2.70 


1992 May 07 




2.2 C 


SBS 1217+499 


AR 


123549.47+591027.0 


2.82 


2006 Jul 29 


7 


3.9 


SBS 1233+594 


PC 


123641.46+655442.0 


3.39 


2006 Aug 29 


7 


4.0 




PC 


135044.67+571642.8 


2.91 


2006 Oct 15 


7 


3.0 


SBS 1348+575 


PC 


140747. 22+645419. 9 


3.08 


2006 Sep 16 


7 


3.7 




PC 


142123.98+463317.8 


3.37 


2006 Sep 15 


7 


3.9 




PC 


142656.17+602550.8 


3.19 


1993 Nov 01 




4.1 c 


SBS 1425+606 


AR, Reimers et al. (1995) 


143835.95+431459.2 


4.61 


2006 Sep 16 


7 


3.7 




PC 


144542.75+490248.9 


3.88 


1993 Jul 12 




5.7 C 




AR 


152156.48+520238.4 


2.19 


2006 Jul 16 


7 


4.1 




PC 


152553.89+513649.1 


2.88 


2001 Dec 08 




24.4 


CSO 755 


AX, BAL, Shemmer et al. (2005b) 


161434.67+470420.0 


1.86 


2003 Dec 20 


4 


2.6 C 


RX J1614.5+4704 


AC, Bade et al. (1995) 


162116.92-004250.8 


3.70 


2001 Sep 05 


2 


1.6 




AC, Bechtold et al. (2003) 


170100.62+641209.0 


2.74 


2000 Oct 31 


1 


39.4 C 


HS 1700+6416 


AC, Reimers et al. (1995) 


173352.22+540030.5 


3.43 


2006 May 28 


7 


3.7 




PC, RLQ 


212329.46-005052.9 


2.26 


2006 Mar 30 


7 


3.9 




PC 


231324.45+003444.5 


2.08 


2005 Sep 23 


6 


1.1 


Hazard 2310+0018 


AC, BAL 


APM 08279+5255 
HS 1603+3820 


3.91 
2.51 


2002 Feb 24 

2002 Nov 29 


3 
4 


H8.8 
8.3 




AC, BAL, Lensed, Chartas et al. (2002) 
AC 



The Chandra exposure lime has been corrected for deleclor dead lime. 

PC — pointed Chandra observation; AC — archival Chandra observation; AR — archival ROSAT observation; AX — archival XMM-Newton observation. For sources with archival X-ray data, we list the paper where the data were originally published, when 
possible. 

C The mean effective exposure time for the aperture used for sources at large off-axis angles. 



threshold of 10 . The probability of spurious detections is 
very low, given our a priori knowledge of the locations of our 
sources and the subarcsecond on-axis resolution of Chandra. 
All Cycle 7 targets were detected by Chandra. 

We report in Table 3 the counts detected in the ultrasoft 
band (0.3-0.5 keV), the soft band (0.5-2 keV), the hard band 
(2-8 keV), and the full band (0.5-8 keV). Manual aperture 
photometry with a 3"-radius aperture was used to derive the 
counts. Also in Table 3 we give the band ratio (the hard- 
band counts divided by the soft-band counts) and the effective 
power-law photon index, Y, assuming an X-ray photon spec- 
trum of the form Ne°^ E~ Y across the full band. This photon 
index was calculated from the band ratio using the Chandra 
PIMMS 6 V3.6A tool; we used the Cycle 7 instrument response 
in PIMMS which accounts for the time-dependent quantum- 
efficiency decay of ACIS at low energies (caused by a thin 
layer of molecular buildup on the ACIS filter). 

We examined the data for the presence of extended X-ray 
emission (e.g., due to gravitational lensing or jets) by com- 
paring the radial profiles of our sources with their expected, 
normalized point-spread-functions (PSFs). None of the Cy- 
cle 7 targets showed extended X-ray emission (but see §2. 1.3). 
We also did not find extended X-ray emission for any of the 
non-gravitationally lensed Chandra archival quasars. 

6 http://cxc.harvard.edu/toolkit/pimms.jsp. 



We checked for an excess of X-ray companions near our 
quasars by searching projected circular regions of w 300 kpc 
in radius centered on each source. We found the number of 
companions within those fields to be consistent with expecta- 
tions from the cumulative number counts from X-ray surveys 
(e.g., Bauer et al. 2004). 

Rapid variability within the observations of our sources was 
searched for by applying Kolmogorov-Smirnov tests to the 
photon arrival times of quasars with > 10 counts. No variabil- 
ity was detected, which is not unexpected given the short ex- 
posure times of our observations ( < 30 min in the rest frame). 

2.3. Archival X-ray Data 

Archival data exist and were utilized for 1 1 of the sources in 
our core sample; in Table 1 we give references to past X-ray 
studies when available. The six quasars SDSS J0145— 0945 
(PI: G. P. Garmire), SDSS J0813+2545 (PI: R. S. Prid- 
dey), SDSS J1614+4704(PI: G. Fossati), SDSS J1621-0042 
(PI: J. Bechtold), SDSS J1701+6412 (PI: L. P. Van Spey- 
broeck), and SDSS J2313+0034 (PI: S. F. Anderson) were 
observed previously with Chandra. With the exception of 
SDSS J2313+0034, all of these sources have > 10 counts 
and, aside from SDSS J1614+4704 and SDSS J1701+6412, 
have been observed on-axis. For SDSS J1614+4704 and 
SDSS J1701+6412, we used apertures with radii of 30" and 
8", respectively, to account for PSF broadening at large off- 
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TABLE 2 

X-ray Properties of Complementary z>4 Sources withM, < -29 



Name 


z 


Total X-ray 
Counts 






Reference 15 


PSS 0133+0400 


4.15 


36 


20.8 


-1.64 


1 


PSS 0134+3307 


4.53 


16 


11.5 


— 1.68 


1 


PSS 0209+0517 


4.14 


22 


15.2 


-1.75 


1 


PSS 0248+1802 


4.43 


14 


21.2 


-1.64 


1 


PSS 0955+5940 


4.34 


10 


6.0 


-1.81 


1 


PSS 0957+3308 


4.20 


17 


11.3 


— 1.74 


1 


PSS 1057+4555 


4.12 


24 


25.3 


— 1.70 


1 


PSS 1347+4956 


4.51 


30 


17.5 


— 1.78 


1 


BR 0241—0146 


4.06 


12 


4.3 


— 1.87 


1 


Tin r\^f\c 1 n ct 


4.73 


3 


2. 1 


— 1.94 


[ 


RR IYV1 1 1799 c 

DJ\ UJ 11 — 1 / ZZ, 


4 00 


7 


4 1 


1 0^ 

— 1 .7J 




BR 0331-1622 


4.36 


14 


8.5 


— 1 .86 


2 


BR 0353-3820 


4.55 


55 


44.7 


— 1 .54 


2 


BR 0418-5723 


4.46 


7 


7.0 


-1.87 


2 


BR 0424-2209 


4.32 


11 


8.6 


-1.81 


2 


PSS 0747+4434 


4.43 


5 


5.7 


-1.81 


2 


PSS 1058+1245 


4.33 


5 


2.4 


-2.07 


2 


BR 1117-1329 


3.96 


2 


1.8 


-2.03 


2 


PSS 1646+5514 


4.04 


5 


5.0 


-2.00 


2 


BR 2213-6729 


4.47 


20 


14.9 


-1.64 


2 


PSS 2344+0342 


4.24 




<2.6 


< -1.98 


2 


Q 0000-263 


4.10 


1229 d 


12.6 


-1.70 


3 


PSS 0926+3055 


4.19 


1156 d 


39.0 


-1.76 


3 


PSS 1326+0743 


4.09 


963 d 


27.9 


-1.76 


3 


BR 2237-0607 


4.56 


306 d 


8.5 


-1.74 


3 



Note. — Quasars are first soi led by reference, then further soi led by right ascension. 

Galactic ab sorption -corrected ilux in the observed 0.5-2 keV band in units of 10 — " erg cm - ^ s — * 

(1) Vignali et al. 2003; (2) Vignali et al. 2005; (3) Shemmer et al. 2005a. 

This source appears in Vignali et al. (2003) under its older designation BR 0308— 1 734. 

Since these sources were observed by XMM-Newton, we quote the total counts from the pn detector. 



axis angles. We averaged the exposure maps over the same 
respective apertures when calculating the count rates. For 
SDSS J0145— 0945, which is a gravitationally lensed quasar, 
we used an aperture with a radius of 4" in order to enclose 
the counts from both images. Although SDSS J0813+2545 
is also gravitationally lensed, the angular separation between 
the images is small enough for a 3"-radius aperture to suffice. 
Data reduction for these sources was carried out in a simi- 
lar manner to that described in §2.2, with the proper correc- 
tion applied to each source for the time-dependent quantum- 
efficiency decay of ACIS at low energies. 

Archival ROSAT data were used for three of the quasars: 
SDSS J1219+4940 (PI: J. Liebert), SDSS J1426+6025 
(PI: D. Reimers), and SDSS J1445+4902 (PI: U. Herb- 
stmeier). These objects were detected by the ROSAT 
PSPC instrument, with 15.5, 8.9, and 28.1 counts, respec- 
tively, in the 0.5-2.0 keV band. SDSS J1219+4940 and 
SDSS J1445+4902 have off-axis angles of 20.6' and 16.1', 
while SDSS J 1426+6025 was observed on-axis. The counts 
were derived for these quasars using aperture sizes of 45" for 
SDSS J1426+6025, 140" for SDSS J1219+4940, and 100" 
for SDSS J1445+4902; background was estimated by placing 
circular apertures on regions where no other X-ray sources 
were present. Using PIMMS to extrapolate our measured 
0.5-2.0 keV count rate for SDSS J1426+6025 to the full 
ROSAT energy range (« 0.1-2.4 keV), we found our value 
was consistent with that previously published by Reimers 
et al. (1995). 

The remaining two archival sources in our sample, 
SDSS Jl 110+4831 and SDSS J1525+5136, were previously 
observed with XMM-Newton. SDSS Jl 110+4831 has its 



X-ray properties described in detail in Page et al. (2005). 
We have re-analyzed the XMM-Newton data using stan- 
dard XMM-Newton Science Analysis System V6.5.0 tasks 
and found results consistent with those presented in Page 
et al. (2005); the X-ray values we quote in Tables 3 and 
4 are taken from our analysis. Similarly, the X-ray data 
for SDSS J1525+5136 have appeared in Page et al. (2005), 
and this quasar has been studied in detail in Shemmer et al. 
(2005b); we use the data from Shemmer et al. (2005b) in our 
analysis. 

3. X-RAY, OPTICAL, AND RADIO PROPERTIES OF 
THE CORE SAMPLE 

In Table|4] (placed at the end of the paper) we list the main 
X-ray, optical, and radio properties of our core sample: 

Column (1 ). — The SDSS J2000.0 quasar coordinates, accu- 
rate to — 0.1". 

Column (2). — Galactic column density in units of 10 20 
cm~ 2 , calculated using COLDEN 7 with the data from Stark 
etal. (1992). 

Column (3). — The monochromatic AB magnitude at a rest- 
frame wavelength of 1450 A (AB^q = — 2.51og / i450 a — 
48.6; Oke & Gunn 1983). The AB1450 magnitudes were calcu- 
lated from the spectra after applying corrections due to Galac- 
tic extinction and fiber light-loss. The fiber light-loss correc- 
tion was calculated as the average difference between the syn- 
thetic g, r, and i magnitudes (i.e., the integrated flux across 
each respective bandpass in the SDSS spectrum) and the pho- 
tometric g, r, and i magnitudes, assuming no flux variation be- 

7 See http://cxc.harvard.edu/toolkit/colden.jsp 
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TABLE 3 

X-ray Counts, Band Ratios, and Effective Photon Indices of the Core Sample 



X-ray Counts" 



Object (SDSS J) 


0.3-0.5 keV 


0.5-2.0 keV 


2.0-8.0 keV 


0.5-8.0 keV 


Band Ratio b 


r b 


U121 JO.U4+144o2i.y 


3 n+2.9 
3 -°-l.6 


TO Q+7.3 

38.8_ 62 


e. a+3.8 
6 - 9 -2.6 


A C 0+7.8 

45 - 8 -6.7 


n 1 e+0.10 


T t+0.5 
z - ;> -0.4 


014516.59—09451 /.J 


c-7 n+8.7 
57 - 9 -7.6 


n+24.8 

566.8I 238 


171 3 + 12.1 

121.3T.n_,, 


/no 1 +27.3 

688.1T 26 . 2 


n 11 +0.02 
°- 21 -0.02 


1+0.1 

2 -!-o.i 


020950. / 1 —000506.4 


< 6.4 


, n n+4.6 
12.0 + 3 4 


c 0+3.6 
3.»_2. 4 


1-7 0+5.3 
I7.8T4.2 


n ,10+0.36 

0.49T .24 


1 3+O.6 
L3 -0.5 


U / JJU2. J l+2Djy 1 1.4 


< 0.4 


33 e+6.9 
-"■°-5.8 


Q+3 8 
6 - 9 -2.6 


An 7+ 7 4 


n 90+0+2 

U - 2U -0.08 


7 7+0.5 
z -0.4 


0/5054.64+425219.2 


3 o+3.2 


45.9+™ 


i n n+4.4 
1U.V_3 3 


56. 8+75 


n T/i+0.10 
°- 24 -0.08 


2 0+° 4 
z - u -0.3 


080342.04+302254.6 


< 6.4 


57.7+^ 


M- 5 3 6 


-+9 5 

71 - 2 -s:4 


°- 23 -ao7 


-» n+0 3 

2 -°-a3 


OfSl JJ1.2O+2545UJ.0 


n+9.2 
66.0_ 81 


A1A 1+21.9 
tjt.i_ 208 


1 Q+13.7 
159.9_ 126 


cqa n+25.4 
594.0_24. 4 


n 3-7+0.04 

°- 37 -0.03 


1 ft+01 


AO/1 /1A1 flC i ficni en n 

084401.95+050357.9 


< 3.0 


1 £ n+5 .1 

16.0_ 40 


-7 0+3.9 
'■ S -2.8 


0+6.O 

23.8T4.9 


n ^o+0.29 
°- 49 -0.21 


1 /i +0.5 


090033 .49+42 1 546.8 


o n+4 
8 -°-28 


O-l n+10 1 

82.OI90 


26.8+j 2 


108.8+1^ 


n 30+0.09 

U - :ii -0.07 


1 t+O 2 

!- 7 -a2 


Oy 4202 . 04 + 042 244 . 5 


£ n+3.6 


3/L Q+7.0 
34.9_ 59 


1 1 7+4J 
11.7_ 34 


Af, 

46 - 6 -6.8 


n 3A+015 
0.34_ 011 


1 7+0-4 


095014.05+580136.5 


< 3.0 


20.8+4J 


i -7+3.5 


26.6+f? 


n To+0.19 
0.2S_ 13 


1 0+O.6 
L8 -0.5 


100129.64+545438.0 


c n+3 4 

5.0+^ 


c . n +8 5 
55.91° j 


10.7+^ 


66.6+^f 


n 1 r»+0 08 

°- 19 -ao6 


1+03 

2J -a3 


A (\A A A>1 1 o i .4 O nm n 1 

101447. 18+430030. 1 


3 -°-i:7 


26.0^5^ 


^ o+3 7 

6 - 8 -2:<5 


32.8+f;* 


n o/;+0 16 

°- 26 -o:!i 


1 n+0 5 

L9 -a4 


1 10610.73+640009.6 


1 -7 n+5 2 
17.014! 


nn A+ll'.O 
"■ 9 -10.0 


23.814:, 


i + 12 2 

123.7+m 


n a +0.07 
°- 24 -0.05 


1 n+0. 2 

L9 Ia2 


11 1038.64+483 115.7 C 












n+0 I 


121930.78+494052.3 d 




15.5+^,' 








1 23549.4 /+59 1027.0 


6.0+ 24 


t i n+6.6 
3L0 -5.6 


<r -7+3.7 


"+7 2 

37 - 7 -6:i 


n ot+0.13 
°- 22 -0.09 


7 0+0.5 

2 -°-0.4 


123641. 46+655442.0 


<4.8 


1 n t+5.5 

19.2+44 


<6.4 


t+5.6 

20.7+45 


<0.33 


> 1.7 


135044.67+571642.8 


<4.8 


i n + 2 fi 

1.91 1 3 


<3.0 


i.7Tf:5 


<1.56 


>2.3 


140/47.22+645419.9 


£ n+3. 6 
6.0+ 2 4 


ia q+6.9 
33.8_ 58 


1 -» 0+4.7 
I2.8T35 


-7+7.9 
46 - 7 -6.8 


n ao+- 16 
O.Jo_„ , 2 


1 <+0.3 
1 - 5 -0.3 


142123.98+463317.8 


<3.0 


t n+3 2 
3 - 9 -L9 


<3.0 


T -7 + 3 1 

3.7+f ; , 


<0.77 


>0.9 


142656. lo+602550.y 




S Q+ 4 + 








OC i Alt' AKCl 1 

143835. 95+43 1459.2 


<o.4 


n s+3.9 
'•°-2.8 


<6.4 


a £+ 4 - 2 


<0.82 


^00 
>0.9 


1 44^4? 764-490248 9 d 




28 l+ 6 4 
zs.i_ 5 3 








152156.48+520238.4 


< 3.0 


<4.8 


1.8+2f 


2 7+ 2 <) 


> 0.38 


< 1.5 


152553.89+513649.1° 








1-81! 


161434.67+470420.0 


5.0«; 4 


126.0+J" 


48.8+^ 


177.8+;^ 


39+ 07 


1 5+ - 2 


162116.92-004250.8 


<4.8 


19 9+ 55 


6 n+3.8 


26.8 + f| 


35+ ' 21 
"■^ -0.15 


1 7+ - 5 
11 '-0.4 


170100.62+641209.0 




273.2+|^j 


74 4 +9 7 


347 6+ 1 ' 9 ' 7 


27+° 04 
u - z/ -0.04 


1 O+0.1 

'• y -o.i 


173352.22+540030.5 


4 0+ 32 
4.U_j 9 


32.8 + |f 


6 Q+3.8 


39 7+ 7 4 
^•'-6.3 


n 2] +0.12 
uzl -0.09 


7 , +0.5 
z - 1 -0.4 


212329.46-005052.9 


<4.8 




<4.8 


23 7+ 5 ' 9 


<0.22 


> 2.1 


231324.45+003444.5 


<3.0 


<4.8 


1.9+f« 


2 q+2.9 


>0.40 


<1.5 


APM 08279+5255 


66.81** 


3967.2+|*^ 


1617.4+-- 


5584.6+^ 


u - 41 -o.oi 


1 3+0.02 
1 - J -ao3 


HS 1603+3820 


8.7i£ 


93.6+ 9 ° 7 7 ' 


23.21=1 


n67 +11.8 
11D -'-10.8 


25+° 07 

u - z;l _0.06 


1 7+ - 2 



a Using Poisson statistics, the errors on the X-ray counts were calculated according to Tables 1 and 2 of Gehrels (1986) and correspond to the 1 O" level. Upper limits on X-ray counts 
were computed according to Kraft et al. (1991) and are at the 95% confidence level. Upper limits of 3.0, 4.8, and 6.4 correspond to finding 0, 1, and 2 X-ray counts within an extraction 
region of I" centered on the optical position of the quasar, considering the background negligible within that region. 

k We define the band ratio as the hard-band counts divided by the soft-band counts. Errors for the band ratio and power-law photon index were calculated at the 1 o level following the 
"numerical method" decribed in §1.7.3 of Lyons (1991). When the number of counts is small, this method avoids the failure of the standard approximate-variance formula. The photon 
indices have been obtained by applying the correction required to account for the quantum-efficiency decay of ACIS at low energies. Note that because the sources in this sample 
have been observed in different Chandra cycles, the quoted band ratios cannot be directly compared with each other, due to the different rest-frame energy bands as well as the time 
dependence of the ACIS quantum-efficiency. 

c The X-ray properties (i.e., T and the soft-band flux) of SDSSJ 1 110+4831 and SDSS J1525+5136 have been calculated from XMM-Newton observations, which observes in a different 
energy band than Chandra. Therefore we do not have values to quote for the X-ray counts. 

d The X-ray properties of these objects were determined from soft-band count measurements taken from archival ROSAT data. Therefore we do not quote values for the X-ray counts 
in the other bands. 



tween the spectroscopic and photometric epochs. The ABuso 
magnitudes for the two quasars that were missed by the SDSS 
selection criteria, APM 08279+5255 and HS 1603+3820, 
were calculated from the photometric i magnitudes. 

Column (4). — The absolute /-band magnitude, taken from 
the SDSS DR3 quasar catalog; for APM 08279+5255 and 
HS 1603+3820, we calculated this value from the photomet- 
ric /-magnitude, correcting for Galactic extinction. 
Columns (5) and (6). — The flux density and luminosity at a 
rest-frame wavelength of 2500 A computed from the magni- 
tude in Column 3, assuming a UV-optical power-law slope of 
a = -0.5 (Vanden Berk et al. 2001), where F v °< v a . 
Columns (7) and (8). — The count rate in the observed- 



frame 0.5-2.0 keV band and the corresponding flux, corrected 
for Galactic absorption and the quantum-efficiency decay of 
Chandra ACIS at low energy. The fluxes have been calcu- 
lated using PIMMS, assuming a power-law model with T = 
2.0, which is a typical photon index for luminous AGNs (e.g., 
Reeves & Turner 2000; Page et al. 2005; Piconcelli et al. 
2005; Shemmer et al. 2005a; Vignali et al. 2005; see §5.1 
for direct justification). 

Columns (9) and (10). — The flux density and luminosity at 
a rest-frame energy of 2 keV, computed assuming F = 2.0 and 
corrected for the quantum-efficiency decay of Chandra ACIS 
at low energy. 

Column (11). — The luminosity in the rest-frame 2— lOkeV 
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band. 

Column (12). 
denned as: 



The X-ray-to-optical power-law slope, Oo X , 



log(/ 2 kev//; 



2500 A J 



log(V 2 kev/V; 



: 0.3838 log(/ 2k ev// 25 oo a) (1) 



2500 i 



where / 2 kev and / 2 5 00 A 816 me monochromatic flux densities 
at rest-frame 2 ke V and 2500 A, respectively. 

Column (13). — The difference between the measured a ox 
(from Column 12) and the predicted a ox (quoted as Aa ox ), 
given the UV luminosity from Column 6, based on the es- 
tablished 0! ox -L 2 5ooA relation [given as equation (2) of S06]. 
The statistical significance of this difference is also given in 
units of <7, where a — 0.146 for 31 <log(L 2500 a)< 32, and 
a = 0.131 for 32 <log(L 2500A )< 33 (see Table 5 of S06). 

Column (14). — The radio-loudness parameter R (see §2.1.3). 
The /4400 A va l ues were calculated by extrapolating from the 
magnitudes in column 3, assuming a UV-optical power-law 
slope of a = —0.5. The /sghz values were calculated us- 
ing data from the FIRST (Faint Images of the Radio Sky at 
Twenty cm; Becker et al. 1995) and NVSS (NRAO VLA Sky 
Survey; Condon et al. 1998) catalogs. The flux density at a 
rest-frame frequency of 5 GHz was calculated from the flux 
density at an observed-frame frequency of 1.4 GHz assuming 
a radio power-law slope of a = —0.8. Seven of our sources 
have FIRST radio detections; two of them are RLQs. Upper 
limits were placed on 20 of our sources at the 3 (7 level, given 
our a priori knowledge of the positions of all of our sources. 
Upper limits on the five sources not covered by the FIRST 
survey were placed using the upper limit on NVSS detection 
(« 2.5 mJy). 

For the luminosity values quoted in Columns 4, 6, 10, and 
11, no lensing corrections have been made (i.e., the fluxes 
have not been de-amplified when determining these values). 

4. OPTICAL SPECTRA AND NOTES ON INDIVIDUAL 
OBJECTS 

In Figure|2]we present the optical spectra of the 32 SDSS 
quasars in our sample. Below we comment on sources with 
particularly interesting optical and/or X-ray properties. 

SDSS J014516.59-09451 7.3 (z = 2.73).— This object was 
first reported as a gravitationally-lensed system by Surdej 
et al. (1987, 1988), and it consists of two images separated by 
A6 ss 2". The optical flux ratio between the two images was 
found to be w 7 (Kassiola & Kovner 1992). Using the Chan- 
dra observations taken in 2004 we have calculated an X-ray 
flux ratio of 12. 4^7 in the full band (see Figure[3]i; this dis- 
crepancy between flux ratios is not wholly unexpected, given 
(1) the amplitude of flux variations in the X-ray band is typi- 
cally greater than that in the optical, (2) the likelihood of in- 
trinsic variability of the quasar between the different epochs, 
and (3) possible microlensing. The optical positions measured 
from SDSS astrometry agree with the positions of the X-ray 
centroids of both components to within the expected astromet- 
ric accuracy of Chandra (w 0.5"). We found a ox = —1.51 for 
this quasar, which differs from the predicted value from S06 
(using a lensing-corrected luminosity) by 1.3 (7. 

SDSS J081331.28+254503.0(z = 1.51).— This object is a 
gravitationally-lensed system consisting of four images dis- 
covered by Reimers et al. (2002). The angular separation be- 
tween the two brightest images is AO = 0.25", too small to 



be resolved by the angular resolution of Chandra (w 0.5"). 
We measure a ox —1.58 for this quasar, which differs from the 
predicted value of S06 (using a lensing-corrected luminosity) 
by 0.4 a. 

SDSS J084401. 95+ 050357.9 (z = 3.35).— This is a BAL 
RLQ, with a relatively mild radio-loudness parameter of R = 
18.9 and a rest-frame equivalent width (EW) of « 30 A for 
the C IV absorption trough. While the additional jet-linked 
X-ray emission generally increases the X-ray-to-optical flux 
ratio, the intrinsic absorption present in BAL quasars tends 
to reduce that quantity; it is possibly because of these com- 
peting effects that we measure oc ox = — 1 .72 for this source, 
which differs from the predicted value for a non-BAL RQQ 
with log(L 2500 A ) = 31.90 at a level of only 0.07 a. 

SDSS J 135044.67+ 57 1642.8 (z = 2.91).— This quasar is 
somewhat X-ray weak, with Oo X = —2.14; this value differs 
from the predicted value by 3.02 a. The SDSS spectrum 
of this object appears in Figure |2j and it shows no obvious 
UV absorption. 

SDSS J152156.48+520238.4 (z = 2.19 s ). — This excep- 
tionally luminous quasar is the third most optically luminous 
object in S05. Manual aperture photometry for this quasar 
measured only 3 counts: 2 in the hard band and 1 in the soft 
band. X-ray fluxes and other properties were calculated from 
the full-band count rate using PIMMS since there were not 
enough counts in the soft band for a detection. This quasar 
is anomalously X-ray weak, with log(vL v )2 kev = 43.74 and 



a steep a 0X = — 2.447IH5, which is inconsistent with the pre- 
dicted value at a level of ss 5 (7. The SDSS spectrum of this 
quasar appears in Figure [4] The Lya line is completely ab- 
sorbed by several narrow absorption-line (NAL) systems, and 
the high-ionization emission lines are blueshifted relative to 
the quasar's redshift, even when allowing for the revised red- 
shift in Footnote 8. The strong observed UV absorption and 
the hint of a hard X-ray spectral shape suggest that absorb- 
ing material along the line of sight is likely responsible for 
the X-ray weakness of this quasar (e.g., Brandt, Laor, & Wills 
2000; Gallagher et al. 2001). Therefore, we will exclude this 
quasar from the statistical analyses below, since our main in- 
terest there is in the intrinsic X-ray emission properties of 
quasars. 

SDSS J170100.62+ 641209.0 (z = 2.74).— This quasar is 
the most optically luminous in the DR3 catalog and has an 
a ox —1.91, which differs from the predicted value by 0.7 <7. 
After binning the data from the X-ray spectrum into 10 full- 
band counts per bin, we fitted the spectrum below 2 keV to 
a power-law model with Galactic absorption and extrapolated 
the fit to higher energies; we found two consecutive bins that 
fall > 3 a below the model near ~ 3 keV in the observed frame 
(~ 11 keV in the rest frame), signs of a possible absorption 
feature. A spectrum taken with XMM-Newton (PI: F. Jansen) 
shows no such feature, although it is not ruled out within the 
uncertainties on the data points (the XMM-Newton spectrum 
was particularly noisy due to background flaring). The binned 
Chandra spectrum of this quasar appears as part of Figure [5] 
Further observations are necessary to test the reality of this 
feature. 



5. DATA ANALYSIS AND RESULTS 



Although the SDSS quotes a redshift of 2.21, upon examination of the 
spectrum we have measured a redshift of 2.19 based upon the Mg II A2798 
line; it is this latter value which we adopt throughout this paper. 
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FLG. 2. — SDSS spectra for our sample of 32 SDSS quasars. Prominent emission lines are marked in the top panel of each column. The spectral resolution is 
~ 2000. The four columns have been sorted into redshift bins in the following way: the first contains sources from z, ~ 1.5-2.5, the second and third columns 
from z ~ 2.5-3.5, and the fourth column from z ~ 3.5—4.5. Within each column the spectra are sorted by right ascension. SDSS J0844+0503, SDSS J1525+5136 
and SDSS J23 13+0034 are BAL quasars. Note the C IV absorption for SDSS J0844+0503, which is also mildly radio-loud, and the interesting spectrum for 
SDSS J1521+5202, which is shown in greater detail in Figure^ 



5.1. X-ray Spectral Properties 
5.1.1. Individual-Object Analyses 

We have investigated the X-ray spectra of seven of the 
eight sources with Chandra observations that have > 100 full- 
band counts (APM 08279+5255 has already had its com- 
plex spectrum studied in detail in Chartas et al. 2002); below 
this threshold there are too few counts for statistically use- 
ful results to be derived from individual spectral fitting. The 
seven spectra were extracted with the CIAO routine PSEX- 
TRACT using circular apertures of 3.0" in radius centered on 
the X-ray centroid of each source, with the exceptions of 
SDSS J0145-0945 (4.0" used to enclose both lensed images), 
SDSS J1614+4704 (30.0" used due to PSF broadening at 
large off-axis angles), and SDSS J1701+6412 (8.0" used due 
to PSF broadening). Background regions were extracted us- 
ing annuli of varying sizes to avoid contamination from other 
X-ray sources. An annulus was not used to extract the back- 
ground for SDSS J1701+6412 due to the quasar's location 



near the edge of the ACIS 13 CCD; instead we chose a nearby 
circular region 20" in radius that was free from other X-ray 
sources. 

We used XSPEC vl 1.3.2 (Arnaud 1996) to fit each spec- 
trum across the full-band energy range (0.5-8.0 keV) with a 
power-law model and a fixed Galactic-absorption component 
(Dickey & Lockman 1990); all fits assumed solar abundances 
(e.g., Anders & Grevesse 1989) and used the WABS absorption 
model in XSPEC. We used the C-statistic (Cash 1979) when 
modelling the unbinned data, since this method is more appro- 
priate when fitting low-count sources than x 2 fitting and still 
remains accurate for higher numbers of counts (e.g., Nousek 
& Shue 1989). All of the errors have been quoted at the 90% 
confidence level considering one parameter to be of interest 
(AC = 2.71; Avni 1976; Cash 1979). Although when using 
the C-statistic there is no value analogous to P(x 2 \ v ) with 
which to perform model testing, we assessed whether each 
model fits the data acceptably by searching for any systematic 
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FIG. 2. — Continued 




SDSS J 1 52 1 56.48+520258.4 z=2. 1 94 



16.80 s 16.70 s 16.60 s 16.50 s 01 h 45 ra 16.40 s 

Right Ascension (2000) 

FIG. 3. — Chandra 0.5-8 keV image of the gravitationally lensed quasar 
SDSS J0145-0945. The image spans ss 7" x 7" on the sky; North is up, and 
East is to the left. The image has been adaptively smoothed at the 2 o~ level. 
The two quasar images are separated by A9 rs 2", and the optical and X-ray 
flux ratios between them are « 7 and ks 12, respectively (note that the optical 
and X-ray flux measurements were taken at different epochs). The positions 
of the two X-ray images agree with those of the optical images (based upon 
SDSS astrometry) to within the expected Chandra positional error. 



residuals. The seven objects with fitting, along with their fit 
parameters and statistics, appear in Table 5. In Figure [5] we 
present their X-ray spectra, binned at a level of 10 counts per 
bin for clearer presentation. Note that in Figure [5] (unlike in 
Table 5) we used % 2 fitting in order to show residuals in units 




1400 1500 
Rest-Frame Wavelength [A] 

FIG. 4.— The SDSS spectrum of SDSS J1521+5202 (see §4). This 
object shows strong Lya absorption, which suggests that absorbing mate- 
rial along the line-of-sight is responsible for its anomalous X-ray weakness 
(0-ox= —2.44). The spectral resolution is rs 2000. 

of a. The values of Y calculated from the band ratios (see Ta- 
ble 3) are consistent with those derived from the best-fit mod- 
els. We also added an intrinsic, redshifted, neutral-absorption 
component to the model, but it did not significantly improve 
any of the fits. 

SDSS J1701+6412 appears to show an absorption feature 
at ~ 3 keV in the observed frame (see §4). We used % 2 fit- 
ting to investigate further the significance of this feature; the 
» 350 full-band counts detected from this object are enough 
for x 2 fitting to be acceptable. When modelled with a power 
law and Galactic absorption, P{% 2 \v) — 0.23 with % 2 = 34.1 
and V = 29. Although this is a statistically acceptable fit, the 
presence of systematic residuals motivated further investiga- 



10 



JUST ET AL. 



TABLE 5 

Individual-Object Spectral-Fitting Results 3 



Object 
(SDSS J) 


X-ray 
Counts 


r b 


& 
Statistic 




yv H c 

(10 22 cm~ 2 ) 


Statistic 


0145-0945 
0813+2545 
0900+4215 
1 106+6400 
1614+4704 
1701+6412 
HS 1603+3820 


686 
591 
108 
122 
181 
352 
116 


? 05+ 012 
Z - U;, -0.12 

1 65+°- 12 

l,y -0.28 

2 ni+o- 28 

z - ul -0.28 
1 74+0.28 
1 - '^-0.22 

1 QI+016 
lyl -0.16 

2 01 +0 29 
ZA -0.28 


151.7 
195.4 
74.5 
73.7 
88.1 
153.3 
63.2 


2.08+°; 2 ° 

1 94+0.30 
'■^ -0.28 

2 oi+ - 37 

ZAn -0.27 

!- 86 -032 

1 91+022 

2 10+ 051 


< 1.00 
<0.16 d 

< 1.56 

< 1.18 

< 1.80 

< 1.47 

< 1.29 


151.4 
195.2 
74.5 
73.7 
88.0 
153.3 
62.8 



a All fits are across observed-frame 0.5-8.0 keV and include appropriate Galactic absorption. 

Without an intrinsic absorption component. 
c With an intrinsic absorption component. Upper limit on the intrinsic column density at the source redshift given in Table 1, quoted at the 90% confidence level. 
^ The tighter constraint on this intrinsic column density is due to the relatively low redshift (z — 1 .5 1) of this object combined with its relatively large number of counts. 
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FIG. 5. — Individual Chandra X-ray spectra and residuals for seven of the 
quasars from our sample which have > 100 counts. The spectra have been 
fitted across the full band (0.5-8 keV). The % residuals are in units of a, and 
the inset in each panel shows a contour map of T versus intrinsic Nu (in units 
of I0 22 cm~ 2 ) at confidence levels corresponding to 68%, 90%, and 99%. 

tion. Performing an F-test, we found the addition of intrinsic 
absorption to the model did not significantly improve the fit. 
However, the addition of an absorption edge at 2.4 keV in 
the observed frame (9.0 keV in the rest frame) significantly 
improved the fit at a confidence level greater than 99.6% 
(Ax 2 = 11.8 for 2 additional fit parameters). Fe XXV has an 
ionization energy of 8.8 keV, which is close to the location of 
the modelled edge. Further observations are required to as- 
sess better the nature of the X-ray spectral complexity in this 
remarkably luminous quasar. 

5. 1 .2. Joint- Spectral Analyses 

The relatively small numbers of counts for many of the 
sources in our sample make it impossible to measure accu- 
rately r and Na on a source-by-source basis, so joint fitting 
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FIG. 5. — Continued 

has been used for different combinations of sources in our 
core sample combined with the complementary z>4 sources 
from S06 (see §2.1.2). This procedure allows measurement 
of the average X-ray properties of sets of objects, including 
those with too few counts for individual spectral fitting. Only 
42 sources that have Chandra observations have been used 
in the joint fitting. Shown in Table 6 are the results of the 
joint fits across the observed 0.5-8.0 keV range, all of which 
exclude gravitationally-lensed (since their deamplified fluxes 
do not meet our luminosity criterion), radio-loud, and BAL 
quasars. Also excluded is SDSS J1521+5202 (see §4), as 
well as sources which had less than 3 full-band counts; these 
include SDSS J1350+5716 (with 2 counts) and two of the 
complementary z>4 sources from S06 (BR 1 1 17—1329, with 
2 counts, and PSS 2344+0342, which was not detected in the 
X-ray band). We extracted the X-ray spectra of these sources 
using PSEXTRACT in a similar manner to that described in 
§5.1.1. The sets of unbinned spectra were fit using XSPEC, 
first with a power-law model and a Galactic absorption com- 
ponent, which was kept fixed during the fit, and then with an 
added intrinsic neutral-absorption component. All fits utilized 
the C-statistic. 

We checked whether our sample was biased by objects with 
a high signal-to-noise ratio (S/N). We split the quasars into 
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two groups, each having more and less than 100 full-band 
counts. It can be seen from Table 6 that the F found for the 
quasars with > 100 counts is consistent with the F found for 
quasars with < 100 counts; we therefore conclude that the 
high S/N spectra do not bias the overall sample (note that the 
total number of counts is approximately split between the two 
groups). The addition of an intrinsic-absorption component 
did not improve any of the fits; i.e., no significant amount of 
intrinsic absorbing material has been detected. Upper limits 
on intrinsic Nh appear in Table 6. 

Previous studies have shown that, in general, F does not 
evolve with redshift for luminous quasars (e.g., Page et al. 
2005; Shemmer et al. 2005a, 2006a; Vignali et al. 2005). To 
investigate this matter further, we have performed joint fit- 
ting on sets of our quasars binned into four integer redshift 
bins ranging from z = 1-5. Our sample of the most optically 
luminous quasars spans the full redshift range where such ob- 
jects are known to exist in the Universe, z ~ 1.5^4.5, and it 
is constrained to a relatively narrow luminosity range [hav- 
ing a mean log(L 2 5oo a) = 32.2; no mean L2500A m anv mte ~ 
ger redshift bin differs from this global mean by more than 
2 <7, with three of the four redshift bins differing by less than 
1 a]. These two properties enable our sample to explore a 
different region of the luminosity-redshift plane than previ- 
ous studies, and they minimize possible confusion between 
redshift-dependent and luminosity-dependent effects. For the 
joint fitting we used the same models described above and 
again found that the high S/N sources (> 100 counts) did not 
bias the best-fit parameters. The best-fit parameters appear in 
Table 6, while a plot of F vs. redshift using these values is 
shown in Figure [6] We found no detectable change in F with 
redshift {% 2 = . 5 for 3 degrees of freedom, P(x 2 \v) —0.91), 
and basic fitting shows that the maximum allowed change in F 
across this redshift range can be no more than ss 5%. We have 
also split the sample into higher and lower optical luminos- 
ity halves and performed joint fitting; within the uncertainties 
the maximum allowed change of F with luminosity can be no 
greater than sa 10%. 

At rest-frame energies below ~ 1 keV, some quasar SEDs 
can deviate from a power-law due to the additional X-ray flux 
provided by the soft X-ray excess (e.g., Porquet et al. 2004). 
This additional soft X-ray emission can bias calculations of F 
toward higher (softer) values. To see if our measurements 
of r are biased in such a way, we redid our joint spectral 
analysis considering only the data above rest-frame 2 keV. 
Since our lower observed-frame energy limit is still 0.5 keV, 
this additional constraint does not affect our calculations at 
z>3. Implementing this additional constraint did not change 
our calculated values of F significantly. As perhaps expected, 
the largest change occurred in the z = 1-2 bin. In this red- 
shift bin, the photon index increased from F = 1 .87^q" }^ to 
T = 1.95 

-0A~i> wmcn i s well within the derived errors (and 
this increase goes opposite to the sense expected if soft X-ray 
excess emission were present). The small difference between 
these two methods shows that our sample is not significantly 
biased by excess soft X-ray emission, so we continue to use 
the results from joint fitting done in the observed-frame, 0.5- 
8.0 keV band. 

We have searched for a narrow, neutral, iron Ka line in 
each of the integer redshift-binned sets of spectra. No lines 
were detected. Upper limits on the rest-frame EWs of any 
such emission appear in Table 7; the rest-frame EWs were 
calculated at the mean redshift in each bin. These upper lim- 
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FIG. 6. — X-ray power-law photon index (F) vs. redshift, binned into 
four integer redshift bins ranging from 7, = 1-5, for our sample and the 
complementary high-luminosity sources at z>4 from S06 (see §2.1.2). The 
mean photon indices were derived from joint spectral fitting across the full 
band (0.5-8 keV). Only sources with Chandra observations are included, 
while gravitationally lensed, radio-loud, and BAL quasars have been ex- 
cluded (as well as SDSS J1521+5202). The vertical error bars show 90% 
confidence bars in each bin, and the number of sources in each bin is shown 
beneath the error bars; three sources with < 3 counts (SDSS 1350+5716, 
BR 1117-1329, and PSS 2344+0342) did not meet our requirements for 
joint fitting and thus were excluded (see §5.1.2). The horizontal error bars 
show the width of each bin, and the data points are marked at the median red- 
shift in each bin. The dashed line shows the best constant-model fit, which 
has a value of 1 .92. 

its range from w 150-500 eV. These are not particularly tight 
constraints, especially when considering the high average lu- 
minosity of our sample. Quasars with higher X-ray luminosi- 
ties tend to show weaker iron Ka emission lines (e.g., Page 
et al. 2004; Bianchi et al. 2007). 

We have also checked for a Compton-reflection continuum 
component at w 10-50 keV in our spectra, which would be 
particularly apparent at high redshifts. No reflection compo- 
nent was found; this is not unexpected given the high lumi- 
nosities and relatively low number of counts for the higher 
redshift sources. 

5.2. X-ray-to-Optical Spectral Energy Distributions 
5.2.1. Basic Sample Properties 

The X-ray- to-optical flux ratio for AGNs has been found to 
decrease at higher optical luminosities, but it does not show 
any clear change with redshift (e.g., Avni & Tananbaum 1986; 
Wilkes et al. 1994; Strateva et al. 2005; S06 and references 
therein; but see Kelly et al. 2007). Using our 34 object core 
sample of highly luminous quasars spanning the widest pos- 
sible redshift range for such objects (z« 1.5-4.5), we fur- 
ther examine the OCo X -L 2500j 4 relationship and provide con- 
straints on dox evolution with redshift. All of the statisti- 
cal analyses presented below have excluded radio-loud, BAL, 
gravitationally-lensed, and weak-line quasars (see below), as 
well as SDSS J1521+5202 (see §4); any group of quasars sat- 
isfying these criteria will be hereafter referred to as "clean". 

Fi sure 171 shows a ox vs. £^500 A f° r our core sample com- 
bined with the full S06 sample (including the complementary 
high -luminosity z>4 quasars) and 14 additional z > 4 quasars 
from Shemmer et al. (2006a), resulting in an X-ray detec- 
tion fraction of 89%. The inclusion of the 14 clean quasars 
from Shemmer et al. (2006a) significantly improves coverage 
at z ~ 5-6; note that the full Shemmer et al. (2006a) sample 
includes four weak emission-line quasars, which we do not in- 
clude in our analyses since the nature of these objects remains 
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TABLE 6 
Joint Spectral-Fitting Results 11 



Sources 


Number of 


Median 


Median 


Total 




& 




7V H C 




Included 


Sources 


Redshift 


Mf 


Counts 


r b 


Statistic 




(10 22 cm- 2 ) 


Statistic 


All RDOs 


42 


4.04 


—29.36 


1872 


1 99+009 


1 149.0 


1 9?+°- 10 


<0.2 


1 149.0 


RQQs, <100 counts 


37 


4.07 


—29.34 


993 


1 Q^+0.10 


695.1 


1 QQ+0.I3 

1 -'-'-o.io 


<0.6 


695.1 


RQQs, >1()0 counts 


5 


2.74 


-29.86 


879 


1 92 +0 Al 
' -0.09 


454.0 


171 -0.08 


<0.8 


454.0 


RQQs, 1 < z < 2 


4 


1.88 


-29.41 


343 


1 87+° ! ? 

— Q. 14 


192.9 


1 91+xH 


<1.2 


192.9 


ROOs 1 < 7 < 3 




2.63 


-29.46 


784 


1 95+0 11 


446.6 


1 95 + !H 7 


<0.8 


446.6 


RQQs, 3 < z < 4 


8 


3.37 


-29.43 


305 


1.90+°;* 

— y. i o 


213.5 


1.86+°" 

— U.18 


<1.0 


213.7 


RQQs, 4 < z < 5 


22 


4.34 


-29.17 


440 


1 Q5+016 
L,y3 -Q.l6 


295.7 


1S7 -0.16 


<0.8 


296.5 


RQQs, <100 counts, 1< z < 2 


3 


1.90 


-29.41 


162 


1 98+ - 25 


103.9 


9 1^+0.48 
iu -0.37 


<2.0 


103.9 


RQQs, <100 counts, 2 < z < 3 


5 


2.82 


-29.39 


194 


i qc+0.23 
1 - yJ -0.23 


156.3 


2 10+ 040 


<2.5 


156.3 


RQQs, <100 counts, 3 < z < 4 


7 


3.38 


-29.41 


197 


, o 7 +0.23 


138.9 


, c, 7 +0.28 
1.8/_ 022 


<2.1 


138.9 


RQQs, <100 counts, 4 < z < 5 


22 


4.34 


-29.17 


440 


, QO+0.16 

'■•"-0.16 


295.7 


i O7+0.T6 
1S/ -0.16 


<0.8 


296.5 


RQQs, High-Luminosity Half 1 


21 


3.70 


-29.57 


1151 


1 Q4+0.09 


681.0 


1 Q1+012 
l,,1 -0.08 


<0.3 


681.2 


RQQs, Low-Luminosity Half 1 


21 


4.03 


-29.14 


721 


1 90+° 13 


468.0 


1 Q7+0.18 
ly '-0.16 


<1.1 


468.0 


RQQs, < 100 counts, High-Lum. e 


18 


4.09 


-29.43 


508 


1 Q3+0.17 


358.9 


1 Q1+0" 
lyl -0.12 


<0.7 


358.9 


RQQs, < 100 counts, Low-Lum. e 


19 


4.06 


-29.02 


485 


1 Q5+0.15 
L,y3 -0.15 


336.3 


1 98+° 23 


<1.3 


336.3 



a Errors on T and upper limits for A^h are quoted at 90% confidence levels. 

Without an intrinsic absorption component. 
c With an intrinsic absorption component. 

The high-luminosity and low-luminosity halves have X-ray luminsities ranging from log(vL v )2 kev — 44.50-45.76 and log(vL v )2 kcv — 44.45^-5.83, respectively. 
e The < 100 count high-luminosity and low-luminosity halves have X-ray luminsities ranging from log(vL v )2 kev — 44.50—45.83 and log(vL v )2 kev — 44.45—45.48, respectively. 



TABLE 7 
Iron Kcc Line Constraints 





Rest-Frame 
EW (eV) a 


Number of 
Counts 


All RQQs, 1< z < 2 


< 490.6 


343 


All RQQs, 2 < z < 3 


< 302.9 


784 


All RQQs, 3 < z < 4 


< 462.0 


305 


All RQQs, 4 < z < 5 b 


< 144.4 


440 


RQQs, < 100 counts, 1< z < 2 


<883.3 


162 


RQQs, < 100 counts, 2 < z < 3 


< 1387.2 


194 


RQQs, < 100 counts, 3 < z < 4 


< 769.7 


197 


RQQs, < 100 counts, 4 < z < 5 b 


< 144.4 


440 



a Calculated at the mean redshift in each bin and quoted at the 90% confidence level. 
b Three of the quasars with < 3 full-band counts were excluded. 

unclear. The best-fit relation from S06: 

Oox = (-0.137 ±0.008)log(L 250oA ) + (2.638 ±0.240) (2) 

is shown as a dotted line; a more detailed analysis of the best 
fit for the correlation between a ox and UV luminosity is given 
in §5.2.2 below. Note that the addition of our core sample of 
26 clean quasars to the full S06 sample increases by a factor 
of ss 2 the number of quasars at the highest luminosities that 
have a ox values. 

In Figure[8]we show histograms of our a ox distribution and 
the distribution of residuals from the S06 best fit (Aa ox ) for 
our 5 1 quasar clean sample. The 26 clean quasars from our 34 
quasar core sample are marked with light shaded histograms, 
while the 25 complementary z>4 sources from S06 are not 
marked. Marked on the Oo X distribution are the measured 
and predicted mean a ox values, as solid and dashed lines, re- 
spectively; the predicted mean a ox value was calculated from 
equation (2). Some standard statistical values for our Oo X and 
AoCqx distributions are shown in Table 8. The mean value of 
a ox = — 1 .80 ± 0.02 for our sample agrees with the predicted 
value from Table 5 of S06 (-1.788) to within 1 a. In the 
AOo X histogram, dotted lines mark the 1 a range from the 



predicted S06 value (see §3, Column 13). As a test, we have 
determined that it would require a reduction of ss 17% to the 
X-ray luminosities of our sources (corresponding to a change 
in a ox of w —0.03) in order for the mean a ox of our sam- 
ple to become inconsistent with the predicted value of S06. 
We have also used the method of Maccacaro et al. (1988) to 
estimate the intrinsic dispersion of the a ox values for our 5 1 
non-BAL RQQs. We find a highly significant intrinsic dis- 
persion of 0.10; the measured dispersion, not correcting for 
measurement errors, is 0. 14. 

As a basic first test for any redshift dependence of a ox , we 
have binned our sample of quasars, including the complemen- 
tary high-luminosity z>4 sources, into integer redshift bins 
from z = 1-5. Recall from §5.1.2 that the mean L 2500 ^ in each 
bin does not differ from the sample mean by more than 2 a, 
with three of the four redshift bins differing by less than 1 a; 
this reduces the effect of the a ox -luminosity correlation when 
looking for any a ox -z correlation. The values for the mean 
0£ ox calculated in each bin, as well as the mean residuals (i.e., 
AoCo X ) from the S06 best fit and the best fit found in this paper 
[see §5.2.2, equation (3)], are plotted against redshift in Fig- 
ure [9] Error bars show the standard error of the mean in each 
bin. No detectable change in oc ox is evident across the full 
redshift range [aconstant fit gives % 2 = 1.3,P(^ 2 |v) =0.74], 
and using basic fitting we have placed a constraint on any such 
change in Oo X to be <6% (corresponding to a change in the 
ratio of L2 kev t° ^2500 A °f l ess tn an a factor of 1 .9). 

5.2.2. Linear-Regression Analysis 

To investigate further correlations between a ox , L 2500 A , 
L2 keV, and z, we have added our core sample of 26 clean 
quasars to the full 333 source sample of S06; we removed 
SDSS J1701+6412 from the S06 sample since it is present 
in both. Also included are the 14 z > 4 clean quasars from 
Shemmer et al. (2006a). The inclusion of our core sample 
of quasars, which lie in a narrow range of high luminosity 
[log(L 2 5ooA) ~ 32.0-32.5] and across a fairly wide range of 
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FIG. 7. — vs. ^2500 A ^ or our core sam Pl e (filled symbols), the full S06 sample, and the Shemmer et al. (2006a) sample (open diamonds); only clean 
quasars have been included. Upper limits are denoted with downward-pointing arrows; all of our core-sample sources have X-ray detections. The solid (dotted) 
line is the best-fit relation found in this paper (S06). The lower panel shows residuals (Aa ox ) from the S06 relation; this equation is shown at the bottom of the 
top panel. The most-luminous object on this plot is PSS 0926+3055, one of the complementary high-luminosity z>4 sources. 



TABLE 8 

ofcx and Aaox Statistical Values 





Number of 




Unweighted 


Measured 


Intrinsic 


1st 




3rd 




Sources 




Mean 


Dispersion 


Dispersion 


Quartile 


Median 


Quartile 


All Objects 


51 


^OX 


-1.797 ±0.019 


0.136 


0.103 


-1.880 


-1.800 


-1.700 






AOo X 


-0.029 ±0.019 






-0.110 


-0.033 


+0.072 


1<Z<2 


4 


^OX 


-1.805 ±0.060 


0.120 


a 


-1.880 


-1.878 


-1.865 








-0.044 ±0.056 






-0.111 


-0.110 


-0.106 


2<z<3 


11 


a ox 


-1.832±0.041 


0.136 


0.093 


-1.942 


-1.841 


-1.758 






AOox 


-0.057 ±0.041 






-0.120 


-0.063 


+0.009 


3<z<4 


11 


^OX 


-1.765 ±0.045 


0.149 


0.125 


-1.880 


-1.713 


-1.678 






AOox 


+0.010 ±0.047 






-0.079 


+0.066 


+0.084 


4<z<5 


25 


Oox 


-1.795 ±0.026 


0.130 


0.070 


-1.871 


-1.796 


-1.709 






Aa ox 


-0.035 + 0.026 






-0.130 


-0.022 


+0.061 



1 Given the low number of sources in this redshift bin, we were unable to accurately determine the intrinsic dispersion. 



redshift (z « 1.5-4.5), allows exploration of a new region of 
the luminosity-redshift plane (see Figure [TOb . Also shown 
in Figure [10] are the additional quasars from Shemmer et al. 
(2006a), which substantially improve coverage at z ~ 5-6. Ul- 
timately, 372 quasars are included in our analysis: 26 from our 
core sample, 332 of the 333 from the full S06 sample, and 14 
from the Shemmer et al. (2006a) sample, increasing the S06 
sample size by « 12%. Note that we do not expect significant 
problematic effects from unidentified BAL quasars at z< 1.5 
in the S06 sample (see §3.3 of Strateva et al. 2005). 

To quantify the correlations found between the optical and 
X-ray properties, we used the Astronomy Survival Analysis 
software package (ASURV rev 1 .2; Isobe et al. 1990; Lavalley 



et al. 1992) to perform linear regressions on the data. ASURV 
treats censored data using the survival-analysis methods pre- 
sented in Feigelson & Nelson (1985) and Isobe et al. (1986). 
We used both the fully parametric EM (estimate and max- 
imize) regression algorithm (Dempster et al. 1977) and the 
semiparametric Buckley-James regression algorithm (Buck- 
ley & James 1979) when performing linear regressions. In 
what follows we report the parameters derived from the EM 
regression, although in all cases the Buckley-James regression 
algorithm agreed within the errors. 

We confirm and strengthen the finding in previous studies 
that a ox decreases with increasing rest-frame UV luminosity. 
Performing linear regressions with ASURV on the combined 
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FIG. 8. — Histograms of (a) a ox and (b) Aa ox for our 51 clean-quasar sam- 
ple, with bin sizes of 0.1 and 0.05, respectively. Our core sample is marked 
(light shaded histogram). The solid (dashed) line in panel a marks the mean 
measured (predicted) a ox of our quasar sample. In panel b, dotted lines mark 
the 1 o~ range (a = ±0.131) for the distribution of AOo X for sources with 
32 < log(L 250Q s) < 33, taken from Table 5 of S06. 

sample of 372 quasars, we found the best-fit relation between 
a ox and L 2500 A to be 

a ox = (-0.140±0.007)log(L 250()A ) + (2.705 ±0.212). (3) 

For comparison, both our best fit as well as the S06 best fit are 
shown in Figure|7]as solid and dashed lines, respectively. 

We also confirm a significant correlation exists between a ox 
and L 2 ke v- The best-fit parameters for this relation are 

«o X = (-0.093 ±0.014)log(L 2 kev) + (0.899 ±0.359). (4) 

Note that the EM and Buckley-James regression algorithms 
are no longer strictly valid when double-censoring is present 
(upper limits exist on both a ox and L 2 keV m the S06 data). 
However, given the high X-ray detection fraction of our com- 
bined sample (89%), we have treated the censored L 2 k e v data 
as though they were detected. 
Studies measuring a relationship between X-ray and UV 

luminosities of the form Lx x have found both j3 ~ 0.7- 
0.8 (e.g., Avni & Tananbaum 1982, 1986; Kriss & Canizares 
1985; Anderson & Margon 1987; Wilkes et al. 1994; Vignali 
et al. 2003; Strateva et al. 2005; S06) and j8 = 1 (La Franca 
et al. 1995). We find the best-fit parameters for the L 2 kev - 
L 25 oo a relation to be 

log(L 2 kev) = (0.636 ±0.018)log(L 250oA ) + (7.055 ±0.553) 

(5) 

while treating L 2 k e v as the dependent variable, and 

log(L 2k ev) = (0.808 ± 0.02 l)log(L 250oA ) + ( 1.847 ±0.694) 

(6) 

while treating L 2 k e v as the independent variable. Using the 
equations given in Table 1 of Isobe et al. (1990), we calculate 
the bisector of the two lines to be 

log(L 2keV ) = (0.709±0.010)log(L 250oA ) + (4.822±0.627). 

(7) 

This result agrees with those previous studies which found j3 
to be inconsistent with unity. 
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FIG. 9. — (a) a ox plotted against redshift, binned into four integer redshift 
bins ranging from z = 1-5. The complementary high-luminosity sources at 
Z>4 from S06 (see §2.1.2) are also included in this analysis. The number of 
quasars in each bin is given below the error bars; these numbers differ from 
those in Figure 6 due to the inclusion of ROSAT and XMM-Newton targets, as 
well as the three sources with < 2 counts (SDSS 1350+5716, BR 1 1 17-1329 
and PSS 2344+0342). The vertical error bars show the standard error of the 
mean in each bin, while the horizontal error bars show the width of each bin; 
the data points are marked at the median redshift in each bin. In panel (b) 
are the residuals from the best-fit relation derived in S06, while in panel (c) 
are residuals from Equation (3) of this paper, i.e. the best-fit relation found 
in our study. The best constant-model fit (giving = —1.80) is shown as a 
dashed line in panel (a). 

Finally, we used ASURV to investigate the relationship be- 
tween Oo X , L 2500 A , and z. We tested three different parametric 
forms of redshift dependence: (1) a dependence on z, (2) a 
dependence on log(l+z), and (3) a dependence on the cosmo- 
logical look-back time, t(z), in units of the present age of the 
Universe. The best-fit parameters for these three relations are 

a ox = (-0.134±0.011)log(L 250oA ) 

- (0.005 ± 0.007 )z + (2.543 ± 0.320) (8) 



a ox = (-0.137 ±0.012)log(L 2500A ) 

- (0.006 ± 0.023)log( 1 + z) + (2.635 ± 0.340) (9) 



«o X = (-0. 143 ± 0.01 l)log(L 2500A ) 

-(0.001 ±0.003)t(z) + (2.824±0.350). (10) 

All three parametrizations have redshift-dependent coeffi- 
cients consistent with zero; note that these equations have the 
same parametric form as equations (8)— (10) of Kelly et al. 
(2007). This finding agrees with previous studies that have 
found no evolution of a ox with redshift (e.g., Strateva et al. 
2005, S06), as well as our results from §5.2.1 (also see Fig- 
ure|9]l. 
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FIG. 10. — (a) ^2500 A an ^ ^ keV vs ' rec lshift for our core sample (26 
filled circles) and the complementary z > 4 sources from S06 (25 filled trian- 
gles). Also shown are the remaining S06 sample (307 open circles) and the 
Shemmer et al. (2006a) sample (14 filled squares). Gravitationally-lensed, 
radio-loud, and BAL quasars have been removed, as well as the extreme out- 
lier SDSS J1521+5202 (see §4). Note the new region of luminosity-redshift 
space being populated by our core sample. 

We have followed the method described in §4.4 of S06 
to compare our results directly with those from earlier stud- 
ies, in particular Avni & Tananbaum (1986), Wilkes et al. 
(1994), and S06. Confidence contours of Aq [the coefficient 
of log(L 2 gooA)] an d [the coefficient of t(z)] were calcu- 
lated using the method outlined in §§3 and 4 of Avni & Tanan- 
baum (1986). Renaming A z as A T , equation (6) of Avni & 
Tananbaum (1986) becomes 

«ox(£ 2 500 A^l x - ra y loud ) = A o[log(L 2 500A) - 30 - 5 ] 

+A t [t(z)-0.5]+A. (11) 

We found best-fit values of [A ,A T .A] = [0.143 ± 
0.011,-0.016 ± 0.041,1.556 ± 0.009] for our sample. 
Figure QT] shows contours of Ao and A T at both the 68% and 
90% confidence levels, along with contours from previous 
studies. Our best-fit value has smaller confidence contours 
than those of previous studies and has an A T consistent with 
zero. 

Some Monte Carlo simulations have suggested that corre- 
lations among a ox , L 250O ^, and L2 keV ma y arise from the 
effects of luminosity dispersion in optically selected, flux- 
limited samples (e.g., Yuan et al. 1998; Tang et al. 2007). 
However, these studies have usually examined the effects 
of luminosity dispersion over a much smaller total range in 
UV luminosity (AlogL 250() ^~ 2.5) than our full sample cov- 
ers (AlogL 250oA ~ 5). In §3.5 of Strateva et al. (2005) the 
authors estimated the dispersions of L 250 o A anc ' ^2 keV (Ouv 
and Ox* when expressed in log units) and, using simulations, 
showed that the dispersions cannot be responsible for the non- 
unity L2keV~ ^2500 A slope they found. They estimated that 
°uv / Ox is not larger than 1 .4 and is plausibly < 1 for their 
sample, and we expect these values to hold for our sample as 
well. The strength of our sample (which builds upon Strateva 
et al. 2005) comes, in part, from the large luminosity range we 
cover, which is much larger than the value of the luminosity 
dispersion in either band. Tang et al. (2007), in their §3, do not 
take full advantage of the large luminosity range and exam- 
ine only a subsample from S06 [and the smaller, high-redshift 
sample from Miyaji et al. (2006)], where the luminosity range 
covered is only slightly larger than the dispersion. Compar- 
ison of the parametric fits calculated by Tang et al. (2007, 
see their Figure 8) with our full sample shows very signif- 
icant disagreement between the data and the fits, especially 
at high luminosities, where our sample of the most-luminous 
quasars helps considerably. Tang et al. (2007), in their §5, 
consider selection effects using Monte Carlo simulations of 



a™ = Ao [log(i uv ) - 30.5] + A, [t'(z) - 05] + A 
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FIG. 1 1 . — Best-fit values and 90% confidence contours for the coefficients 
Ao and A T for our sample of 59 quasars combined with the S06 sample and 
the Shemmer et al. (2006a) sample (cross, dark solid contour), a total of 
372 clean quasars. Also shown are best-fit values and contours for the S06 
sample (cross, light solid contour), the sample of Avni & Tananbaum (1986; 
open circle, dot-dashed contour), and the sample of Wilkes et al. (1994; 
filled circle, dashed contour). Inset: Magnified view of the 68% and 90% 
contours for both our sample (solid contours) and the S06 sample (light solid 
contours). 

a sample more similar to the full S06 sample. They assume 
the true slope of the L2 kev~ ^250oA relation is unity and as- 
sess how much the observed slope, /3, can differ from unity. 
They require unrealistic combinations of Guv and Ox (e.g., 
with Otjv/Ox 1) in order to obtain values of j3 as flat as 
the 0.709 that we measure in equation (7). While our /3 value 
may be biased somewhat downward owing to inevitable se- 
lection effects in the presently available samples, it seems un- 
likely that the true slope of the L2 keV - ^2500 A relation could 
be unity for optically selected quasar samples. 

6. CONCLUSIONS AND SUMMARY 

We have analyzed the X-ray properties of a large and sta- 
tistically representative sample of the most optically luminous 
quasars spanning M, « —29.3 to —30.2 across a redshift range 
of z ~ 1 .5^1.5. Our total quasar sample consists of 59 sources, 
including 32 from the SDSS, two additional quasars that were 
missed by the SDSS selection criteria, and 25 comparably lu- 
minous z>4 quasars. All of these sources have sensitive X-ray 
coverage from either targeted Chandra observations or from 
archival Chandra, ROSAT, or XMM -Newton observations; 58 
of the 59 sources (98%) have X-ray detections. For some of 
our analyses we have included 332 quasars from S06 and 14 
quasars from Shemmer et al. (2006a). Our main results, de- 
rived for radio-quiet, non-BAL quasars that are not gravita- 
tionally lensed (i.e., our "clean" quasars), are the following: 

• The mean X-ray power-law photon index for our sam- 
ple of the most-luminous quasars is T = 1 .92^Jj {jg, con- 
sistent with values found in earlier studies. 

• Any intrinsic absorbing material for the most-luminous 
quasars has been constrained to have a mean column 
density of Nn<2 x 10 21 cm~ 2 , showing that the most- 
luminous quasars typically have little intrinsic X-ray 
absorption. 
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• Using joint spectral fitting, we have found no significant 
change in T with cosmic time over the redshift range 
z » 1.5-4.5. 

• The mean Oox value of the most-luminous quasars is 
Oqx = —1.80 ±0.02; this is in agreement with earlier 
studies and agrees with the predicted value from S06 
[at a mean luminosity of log(L 2 5oo a) — 32.2] to within 
la. 

• In our sample there is no significant change in 0Co X with 
redshift from z « 1 .5-4.5 when binned and compared to 
a constant model, consistent with results found in some 
earlier studies. 

• Combining our sample with that of S06 and Shemmer 
et al. (2006a) results in a clean sample of 372 quasars. 
Using a parametric modelling method on this sample, 
we found a ox is clearly dependent on L 2500 ^, but shows 
no significant dependence on redshift (for three dif- 
ferent parameterizations of redshift dependence). The 
X-ray-to-optical flux ratios of quasars have not signif- 
icantly evolved out to z ~ 6, and in particular have 
not significantly evolved out to z ~ 4.5 for the most- 
luminous quasars. 
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X-ray, Optical, and Radio Properties of the Core Sample 
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a Flux density al resl-frame wavelenglh 2500 A in units of 10~ 2 ^ ergs cm - 2 s — * Hz - 

k Observed count rate computed in the 0.5-2 keV band in units of 10 — ^ counts s — ' . The count rates for off-axis sources have been corrected for vignetting. 
C Galactic absorption-corrected flux in the observed 0.5-2 keV band in units of 10~ 15 erg cm -2 s~ 1 . 
^ Flux density at resl-frame 2 keV in units of 10~ 32 ergs cm -2 s -1 Hz -1 . 

The difference between measured and predicted a ox {Aa ox )- ™d the significance of dial difference ( a), based on the S06 a ox -L v (2500 A) relation. 

f 

Flux density Lit an ohserved-iraine [reqLiency ol 1 .4 Gil/ taken from fhe NVSS survey. All of her 1 .4 Gil/, flux densities are from Lhe FIRST survey. 
^ Actox values for grav nationally lensed objects have been calculated using lensing-corrected luminosities. 



